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A standard and reasonable definition of a symptotic veloc ity term dominance (AVTD) shows that 



Gowdy cosmologies on X R have an AVTD singularity. 



the numerical study by Hern and Stewart ( 5r-qc/9708036 ) confirms previous results that generic 



The wave equations for Gowdy universes on x R can be written as 

P,rr -e^'^Ql = e-2-(p,g, -e^^g,^ ) (1) 



and 



2P,,Q,,= e-^"(Q,ee+2P,eQ,e) (2) 



where P and Q are respectivelyproportional to the amplitudes of the + and x gravitational wave polarizations and 
,a signifies d/da. One usually Q| defines asymptotic velocity term dominance (AVTD) to mean that, as r — > oo, P 
and Q approach a solution to the equations obtained by setting the right hand sides of (|l|) and (||) equal to zero. 
These "AVTD equations" have the exact "AVTD solution" ||] 

P = ln[ae-^'^(l + C^e^^'^)] (3) 

Q= ^J(T$S-^+e (4) 

with the limiting behavior as r — > cxd of 

P = vT, Q^Qo (5) 

where a, Ci ^ are arbitrary functions of 0, u > is a function of 9, and Qo = l/aC + Grubisic and Moncrief 
(GM) have argued |^ that, given the limiting AVTD behavior (||), the spatially dependent term e^^Q,l in (|^) would 
"spoil" the AVTD solution if |ti| > 1. Since ?; > is already required, GM conjecture that, for generic Gowdy models 
to approach the AVTD solution as r — > oo, one requires < v < 1 everywhere in 9 asymptotically except, perhaps, 
at a set of measure zero. If the Gowdy models have an AVTD singularity, one would expect P and Q to approach 
and (4) with < w < 1 (almost) everywhere for sufficiently large r. This behavior has been confirmed in our 
numerical simulations and has, in fact, been confirmed by HS. To obtain the AVTD solution for a given system 

of equations, one sets the spatial derivative terms to zero by hand. However, our claim on the magnitude of the 
spatial derivative terms in an AVTD regime is merely that the dynamics (character of the solution) is that of the 
AVTD solution so that the spatial derivative terms do not influence the dynamics. Using this standard interpretation 
of an AVTD singularity, one sees that HS's discussion of the magnitude of spatial derivative terms is irrelevant. In 
fact, one can use GM's method of consistency to explain HS's Fig. 3. This is easy to see if one defines ttq = e^^Q,r- 
Substitution in the AVTD equations shows that ttq is a constant (in r) in the AVTD regime |]. Assume the AVTD 
limiting form (||) for P and Q and substitute these in (Q). One finds 

P,rr = -e-2[i-"(^)l^[Q^,(0)]2 + e-^-v"{9)T + 7r|(0)e-2-W- (6) 

where ' = d/d9. GM's conjectured restrictions on the magnitude of v mean that either the first or third term will 
dominate as r — > oo depending on whether v is less than or greater than 1/2. This effect, entirely consistent with the 
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presumption of an AVTD singularity in the Gowdy models, allows one to reproduce HS's Fig. 3 by plotting regions in 
the 6-T plane where v is greater than or less than 1/2. Finally, rewriting (||) as 

^Q,r = e-'"(e'^Q,9),e (7) 
shows consistent behavior in an AVTD regime since the right hand side becomes 

g-2[l-.(9)]r[g/,(^) + 2v\e)Q'^(e)T] (8) 

which decays faster or slower than the third term on the right hand side of (^) depending on the sign of w — 1/2. 

HS also criticize previous results |^,||j6j describing the approach to the AVTD solution. We have argued that the 
nonlinear terms in (|l|) act as potentials Vi = c'^^ttq and V2 = e~^''^~^)(3,g to drive the system to the AVTD regime. 
GM's conjectured generic AVTD behavior requires (as t 00) P,t- > and \P,r \ < 1 since P,t v asymptotically. 
As was pointed out elsewhere [^,^, if P,t < 0, Vi becomes important eventually causing a bounce which changes the 
sign of P,T- On the other hand, if P,,- > 1, V2 will become important eventually causing a bounce which changes the 
sign of P,r —1. (If this causes P,r to become negative, there will be another bounce off Vi and so on.) Both potentials 
will become negligible when the AVTD limit is reached. Clearly, non-generic behavior can occur if either ttq or Q,g 
is precisely zero so that the effect of the potential is absent. For these coefficients approximately zero, it takes a long 
time for the respective bounce to occur since the potential is quite flat. As one increases the spatial resolution, the 
simulation grid points are closer to the non-generic points. Thus the potential is flatter at such points than it would 
be at nearby grid points within a coarser resolution. This leads to the effect, as shown in H,^, that it takes longer for 
a higher resolution simulation to appear to be AVTD (with < u < 1) everywhere, exactly as was found by HS. For 
example, Fig.l of |^ displays the maximum value of v on the spatial grid vs r for two spatial resolutions. The coarser 
resolution was run sufficiently far to yield v < 1 everywhere (in constrast to HS's claim that this was not done). The 
finer resolution simulation was only run far enough to indicate the expected resolution dependence. This can also 
be seen by comparing P,t vs at a given r for different spatial resolutions. One finds fewer points with w > 1 in a 
coarser simulation. Careful examination of the waveforms at 9 values where v > 1 shows that, at every such point, 
Q has an extremum in 9 as one expects in our picture and as has been described in 1^,^,^ . The non-generic behavior 
associated with ttq « leads to apparent discontinuities in Q as has been discussed elsewhere . 

In 1^, we show that, at a fixed value of r, there exists a spatial resolution which will resolve any given peak in the 
waveforms but that, at later t, the peak has narrowed so that greater spatial resolution is needed to resolve it. For 
example, a peak near Q,0 = will narrow as t increases because (P,^- —1) changes sign first at 9 values farthest from 
the center of the peak. Thus P will continue to grow at the original v > 1 only at points ever closer to the center. 
This interpretation of our results agrees with the fact that the equations (^ and (||) obtained by assuming the 
AVTD solution as r — *■ 00 are inconsistent if v{9) > 1 [|j unless both Q'q and Qq vanish. However, one gencrically 
expects a set of measure zero at which Qq = but Qq ^ 0. Thus, using the strict definition of AVTD given by 
Isenberg and Moncrief one concludes that our numerical results support the view that generic Gowdy models 
have a singularity that is AVTD everywhere except, perhaps, at a set of measure zero. To be consistent with our 
results, we would expect that an adaptive mesh refinement (AMR) algorithm would easily observe this narrowing of 
the spikes and insert ever more points about the spike center as t increases. However, HS do not find this behavior 
rather finding that the peaks have "softened somewhat" with v < 1 everywhere. Since the Lax-Wendroff scheme is 
known to lead to a decrease in the amplitude of peaks, we contend that HS's AMR code is not correctly modeling the 
peaks at late r. 

Finally, we wish to point out that the Gowdy simultions described in [|[||J|] used the symplectic algorithm described 
in detail in Given sufficient spatial resolution, this algorithm is stable for evolution of the Gowdy equations. (In 
, the resultant data were smoothed prior to graphing in some of the figures. This "averaging" had nothing to do 
with the simulation algorithm.) 

In summary, the AMR algorithm of HS has generated a Gowdy evolution that confirms previous analytic and 
numerical l2|j|,P results except close to the centers of peaks where we believe our algorithm to be more acurate. 

ACKNOWLEDGEMENTS 

B.K.B. and V.M. would like to thank the Albert Einstein Institute at Potsdam for hospitality. This work was 
supported in part by National Science Foundation Grants PHY9507313, PHY9722039, and PHY9503133. 



2 



[1] Isenberg,J., Moncrief, V. (1990) Asymptotic Behavior of the Gravitational Field and the Nature of Singularities in Gowdy 

Spacetimes Ann. Phys. (N. Y.) 199, 84122. 
[2] Berger, B.K., Moncrief, V. (1993): Numerical investigation of cosmological singularities. Phys. Rev. D 48, 4676 
[3] Grubisic, B., Moncrief, V. (1993): Asymptotic behavior of the x R Gowdy space-ti mes. Phys. Rev. D 47 2371 
[4] Hern, S. J., Stewart, J.M. (1997) The Gowdy Cosmologies Revisited. ^r-qc/97080"3i 

[5] B.K. Berger (1996) "Numerical Investigation of Cosmological Sing ularities" in Re lativity and Scientific Computing, edited 

by F.W. Hehl, R.A. Puntigam, H. Ruder (Springer- Verlag, Berlin), [gr-qc/951200^ . 
[6] B.K. Berger (1997) "Numerical Investigation of Singularities" in Proceedings of the 14th Internatial Conference on General 

Relativity and Gravitation, edited by M. Francaviglia, G. Longhi, L. Lusanna, E. Sorace (World Scientific, Singapore), 

gr-qc/9512003|. 



3 



